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Abstract 


Introduction 


A source term model for steads micro jets was 
implemented into a non-proprietary Navier-Stokes 
computer code. OVERFLOW. The source term 
models the mass flow and momentum created by a 
steady blowing micro jet. The model is obtained by 
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jet to the Navier-Stokes equations. 

The model was tested by comparing with data 
from numerical simulations of a single, steads micro 
jet on a flat plate in tsso and three dimensions. 

The source term model predicted the velocity 
distribution ssell compared to the tsvo-dimensional 
plate using a steads mass floss boundary condition, 
ss hich ssas used to simulate a steady micro jet. The 
model ssas also compared to tsso three-dimensional 
flat plate cases using a steads mass floss boundary 
condition to simulate a steads micro jet. The three- 
dimensional comparison included a case ssith a grid 
generated to capture the circular shape of the jet and 
a case without a grid generated for the micro jet. The 
case svithout the jet grid mimics the application of the 
source term. The source term model compared ssell 
ssith both of the three-dimensional cases. 
Comparisons of selocits distribution ssere made 
before and after the jet and Mach and sorticity 
contours ss ere examined. 

The source term model allosvs a researcher to 
quickly investigate different locations of indisidual 
or ses eral steads micro jets. The researcher is able to 
conduct a preliminary ins estigation ssith minimal 
grid generation and computational time. 
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Boundary layer ingesting (BLT) serpentine (S- 
duct) inlets have the potential benefits of reduced 
drag. size, and sseight by eliminating the boundary 
laser diverter and shortening the inlet duct. Tsso 
drasvbacks of BL1 S-duct inlets are increased 
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recosers. Loss-profile vortex generating (VG) vanes 
hase shossn significant improvement in turbofan 
engine-face distortion in the design of compact 
aircraft inlets [1], Vortex generating vanes can 
increase the radar signature of an inlet and can also 
be difficult to maintain. The use of non-intrusive 
vortex generating techniques such as micro jets, are 
being investigated at NASA Langley Research 
Center by Gorton et al. [2] and other institutions. 

Computational fluid dynamics (CFD) is a tool 
that can be used in conjunction ssith experiments to 
design ness inlets. Before CFD can be used, 
researchers have to be sure that CFD analy sis can 
model the complex flows associated ssith the vortices 
created bs a vortex generating des ice such as vanes 
and micro jets. Allan et al. [3] hase shossn a 
comparison betsseen a CFD analysis and 
experimental data of a loss-profile vortex generator 
s ane on a flat plate. Their results indicate that CFD 
analysis can accurately predict various aspects of the 
resulting vortex including size, shape, location, and 
decay. While the svork conducted bs Allan et al. ssas 
significant, the method used to model the sane has 
mans drass backs. In particular, the small size of the 
vane required vers fine grids to model the vane and 
the region immediately surrounding the vane. The 
result ssas an increase in the number of grids required 
to model the entire floss and ssith an increase in grids 
comes an increase in computational time. Modeling 
indis idual micro jets also results in an increase in 
grid size and computational time. As the number of 
vortex generating devices is increased from one to 
possibls tsventy. a more efficient method must be 
used. 


A source term model for an arras of vortex 
generating vanes ssas des eloped by Bender et al. [4] 
and implemented by Waithe [5] into the Reynolds- 
Averaged Navier-Stokes (RANS) CFD code 
OVERFLOW version 1.8s [6], svhich is a non- 
proprietary CFD code. The s ane source term model 
ssas the first of three models to be implemented into 
OVERFLOW. Once complete the user ss ill has e a 


choice of using a model for vanes, steady micro jets, 
and pulsing micro jets. These models will allow the 
user to quickly investigate vortex generating devices. 

This paper describes the implementation of the 
source term model for steady blowing micro jets in 
the Reynolds Averaged Navier-Stokes (RANS) CFD 
computer code OVERFLOW version 1.8s. In 
addition, this report presents the results of CFD 
analysis utilizing the source term model on a two- 
dimensional (2D) flat plate modeling a single steady 
micro jet. The results are compared to a CFD 
analysis where a constant mass flow boundary 
condition is used at the surface to model the micro 
jet. Also, CFD analysis using the source term model 
to model a steady micro jet on a three-dimensional 
(3D) flat plate is compared to two computational 
cases of the same flat plate with a steady mass flow 
boundary condition to model a steady micro jet. The 
first case utilizes a grid generated to capture the 
shape of the micro jet. The second case does not 
utilize a grid to capture the shape of the micro jet, but 
applies the mass flow boundary condition over a 
rectangular region of the existing grid. This mimics 
the application of the source term model. 

Nomenclature 

A cross-sectional area 

A, cell area 

Ajc, jet area 

A lot total cell area of source term application 
BC48 steady mass flow boundary condition 
c calibration constant 

dA differential cross-sectional area 

K cross-stream kinetic energy 

M Mach number 

MFR JCI jet mass flow rate 
Sj mass flow term 

S, momomentum term 

S 5 energy term 

t unit vector normal to jet 

t, component i of vector t 

u component 1 of U 

U-, free stream velocity 

U local velocity vector 

v component 2 of U 

w component 3 of U 

Wjc jet velocity 

x axial direction 

y cross stream direction 

z normal direction 

d angle of attack of jet to local flow 

d angle of incidence of jet to surface 

d density 


d s free stream density 

d stream-wise vorticity 

civ differential of v 

dw differential of w 

dy differential of y 

dz differential of z 


Abbreviations 

2D two-dimensional 

3D three-dimensional 

BLI boundary layer ingestion 

CFD computational fluid dy namics 

RANS Reynolds Averaged Navair-Stokes 

SA Spalart and Allmaras 

VG vortex generator 

Numerical Method 

The basic methodology of the model is to 
introduce mass flow, momentum and energy terms to 
the mass flow, momentum, and energy equations 
representative of the mass flow, momentum, and 
energy that are added by a steady blowing micro jet. 
The terms are given in equations 1 , 2 and 3. 
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In equations 1 to 3, S, is the mass flow term, S, is the 
momentum term, S 5 is the energy term, c is a 
calibration constant. A, is the cell area, A u „ is the total 
cell area, A jct is the area of the jet, d,, is the free 
stream density, U„ is the free stream velocity, MFR jeI 
is the mass flow rate of the jet, and t, is the directional 
component of the unit vector t . 

The source term model was implemented into the 
CFD code OVERFLOW version 1.8s. OVERFLOW 
solves the steady, compressible, Reynolds-Averaged 
Navier-Stokes (RANS) equations using the diagonal 
scheme of Pulliam and Chaussee |7|. The RANS 
equations are solved on structured grids using the 
overset grid framework of Steger et al. [81, which 
allows for complex geometries. The source term 
model is invoked in OVERFLOW by specifying a 
region of cells containing the steady micro jet to be 
modeled as shown in figure 1. In particular, the cells 
that enclose the area of the jet and the second row of 
cells from the surface are selected to add the 
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momentum and mass flow of the jet. The first row is 
not selected because the source term model is not a 
boundary condition and the first row of cells in 
OVERFLOW is reserved for boundary conditions. 
Since the cells are typically very dense in the 
boundary layer, the application on the second row of 
cells is sufficient. 


Top View 




Figure 1 . Schematic of cells selected for application 
of source term model. 

The user also specifies the calibration constant c. jet 
mass flow rate. MFR iet based on the free stream flow 
conditions, the jet area. A, et . the angle of attack 
between the local flow and the micro jet. d , and the 
angle of incidence to the surface. 8. Figure 2 show s a 
schematic w ith the orientation of the jet to a surface. 
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Figure 2. Schematic of jet orientation. 


2D Flat Plate Study 


The source term model was tested by running 
simulations on a 2D flat plate using the source term 
to model a steady micro jet. The model was 
compared to simulations on a 2D flat plate using a 
steady mass flow boundary condition to model a 
steady micro jet. The grid for the flat plate used in 
this study is show n in figure 3 and comprises a single 
block grid with a total of 9,333 grid cells. 



Figure 3. Computational domain forlD flat plate 
study. 

The grid is dimensionless. The jet modeled has an 
area. A |C1 =0.015. The source term model was applied 
with an angle of attack. 8=0° and an angle of 
incidence. 8 = 0°. If the angle of incidence or angle 
of attack w ere not zero, a grid for the jet w ould have 
been generated for use with the steady mass flow 
boundary condition in order to generate the desired 
angles of attack and incidence. The steady mass flow 
boundary condition does not allow the user to specify 
an angle of attack or incidence: they are both zero by 
default. The source term model allows the user to 
investigate varying angles of attack and incidence 


3 





with out having to generate a grid for the jet being 
modeled. All simulations were run at a free stream 
Mach number, M„=0.2, a Reynolds number, 
Re= 1 x 1 0\ and a jet mass flow rate, MFR JC1 =1.0. 

The multi-grid acceleration option of 
OVERFLOW was used to improve the convergence 
of the steady-state solution. All simulations were run 
using the one-equation Spalart and Allmaras (SA) 
model |9|. The SA model was used because of its 
simplicity. 

OVERFLOW was run on a Compaq Alpha 500 
MHz machine for this study. A total of nine runs 
were completed with the source term model to 
investigate the effects of the calibration constant and 
one run with the steady mass flow boundary 
condition. The total run time for each case was 
approximately one minute. 


2D Flat Plate Results 


The results of the source term model simulations 
were compared to simulations using a steady mass 
flow boundary condition. The source term model was 
calibrated before comparisons were made. The 
model was calibrated by plotting the effects of the 
normalized, cross-stream kinetic energy, K, 
downstream of the jet as the calibration constant is 
increased. As the calibration constant is increased, 
the cross-stream kinetic energy downstream of the 
source term should asymptote to a constant value. 
The source term model is correctly calibrated by 
using values of the calibration constant that occur in 
the asymptotic region of the cross-stream kinetic 
energy plot. The normalized, cross-stream kinetic 
energy is given by equation 4. 
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In equation 4, u, v. and vv are the components of the 
local velocity vector, U , d is the density of the flow 
field, and A is the area of the cross-section over 
which the cross-stream kinetic energy is integrated. 
Figure 4 shows the square root of the normalized, 
cross-stream kinetic energy, K 12 , versus the 
calibration constant, c. For comparisons with the 
steady mass flow boundary condition, a calibration 
constant of c=3.8 was used. 


0.038 
0.036 
0.034 
0.032 
N 0.030 
* 0.028 


0.026 

0.024 

0.022 

0.020 


♦ Source Term 
- - BC48 


1.0 1.5 2.0 2.5 3.0 3.5 4.0 

C 


Figure 4. Square root of cross-stream kinetic energy 
vs. calibration constant, c. 


Figure 5 and 6 shows the mach contours of the 
source term and steady mass flow boundary 
condition, respectively. 
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Figure 5. Mach contours using the source term 
model. 
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Figure 6. Mach contours using the steady mass flow 
boundary condition. 
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Figures 7 and 8 show close ups of the jet with Mach 
contours and velocity vectors of the source term 
model and steady mass flow boundary condition, 



respectively. The source term model does a ven 
good job in reproducing the effects of the stead} 
mass flow boundary condition. 



0.13 

0.17 

C.15 
0.13 
0.1 1 
0.03 
0.07 
0.04 
0.02 
C.GC 


Figure 7. Close up of jet with Mach contour using the 
source term model. 



Figure 8. Close up of jet with Mach contours using 
the steady mass flow boundary condition. 


The velocity profile at four different locations was 
examined to evaluate the effect the source term 
model had on the flow field compared to the steady 
mass flow boundary condition. Figures 10-13 show 
the velocity profiles caused by the source term model 
and the stead} mass flow boundary condition. In 
figures 9-12. the steady mass flow boundary 
condition is labeled as BC48. which is its 
OVERFLOW name. At x=0.05 the velocity profiles 
of the source term and steady mass flow boundan 
condition are identical, which indicates the source 
term model and stead} mass flow boundary condition 
produce similar effects upstream of the jet. The 
velocity profiles are only slight!} different at 
locations immediate!} after the jet and very close to 
the surface. As the flow field moves away from the 
surface and away from the jet. the profiles are 
identical. The discrepancy between the profiles close 


to the surface could be attributed to the different 
locations w here the source term and steady mass flow 
boundan condition are applied. In particular, the 
stead} mass flow boundan condition is applied at the 
surface, while the source term is applied one cell 
aw ay from the surface. 
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Figure 9. Velocity distribution at x=0.05. 
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Figure 1 0. Velocity distribution at x=0.15. 
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Figure 11. Velocity distribution at x=0.1 7. 
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Figure 12. Velocity distribution at x=0. 19. 

A representative flow residual history of both the 
source term model and steady mass flow boundary 
condition is shown in figure 13. 
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Figure 13. Flow residuals history. 


The residual history is used to track the convergence 
of the CFD solution. A residual drop of several 
orders of magnitude is expected from the 1/4 grid 
resolved solution to the full grid resolved solution of 
a converged solution. The solution has dropped close 
to three orders of magnitude. 


3D Flat Plate Study 


The source term model was also tested by 
running simulations on a 3D flat plate using the 
source term to model a steady micro jet. The model 
was compared to two simulations on a 3D flat plate 
using a steady mass flow boundary condition to 
model a steady micro jet. One of the simulations 
using the steady mass flow boundary condition used 
additional grids to capture the shape of the micro jet. 
The other simulation using the steady mass flow 


boundary condition did not use additional grids to 
capture the shape of the micro jet. 





Figure 14. Computational domain for3D flat plate 
study. 



Figure 15. Close up of computational domain for 
steady mass flow boundary condition with jet grid. 



Figure 16. Top view, close up of computational 
domain for steady mass flow boundary condition with 
jet grid. 

The boundary condition was applied to rectangular 
grid cells similar to the application of the source term 
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model. The grids used in this study are shown in 
figures 14-16. The grid used for the source term 
model and the steady mass flow boundary condition 
with out the jet grid comprises a single block grid 
w ith a total of 698.139 grid cells. The grid used for 
the steady mass flow boundary condition with the jet 
grid comprises the single block grid used by the 
source term model (698.139 grid cells), a cylindrical 
grid for capturing the shape of the jet (48.461 grid 
cells) and an axis grid (2.783 grid cells) for a total of 
749,383 grid cells. The axis grid is necessary to 
remove the singularity created by the cylindrical grid. 

The grids are dimensionless. The jet modeled 
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was applied with an angle of attack, d = 0° and an 
angle of incidence, b = 0°. All simulations were run 
at a free stream Mach number. M„=0.2. a Reynolds 
number. Re=lxl0\ and a jet mass flow rate. 


MFR jeI =l .0. 

The multi-grid acceleration option of 
OVERFLOW was used to improve the convergence 
of the steady-state solution. All simulations were run 
using the one-equation Spalart and Allmaras (SA) 
model [9]. The SA model was used because of its 
simplicity. 

OVERFLOW was run on a Compaq Alpha 500 
MHz machine for this study. A total of three runs 
were completed with the source term model and the 
steady mass flow boundary condition. The total run 
time for each case was approximately ten minutes. 


3D Flat Plate Results 


The results of the source term model simulations 
were compared to two simulations using a stead} 
mass flow boundary condition to model a stead} 
micro jet. One of the simulations using the steady 
mass flow boundary condition utilized a grid for 
capturing the shape of the jet. w hile the other did not. 

Figures 17-19 show the Mach contours of all 
three cases at the centerline of the jet, y=0. The 
source term model has features of both of the 
simulations with the steady mass flow boundary 
condition. In particular, the jet looks more like the 
simulation w ith the grid generated for the jet. but the 
region immediately follow ing the jet looks more like 
the simulation w ithout the grid generated for the jet. 

Figures 20-23 show the centerline u velocity 
distribution for all cases. The stead} mass flow 
boundary condition is noted as BC48 on the plots. 
The velocity distributions are identical at x=0.985 as 
expected, which indicates very the three methods 
produce similar effects upstream of the jet. 


Immediately after the jet, at x=1.015, the source term 
model has a profile that is similar to the steady mass 
flow boundary condition with the jet grid close to the 
surface (z=0). 
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Figure 17. Mach contours using the source term 
model. 



X 


Figure 18. Mach contours using the steady mass flow 
boundary condition without jet grid. 
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Figure 19. Mach contours using the steady mass flow 
boundary condition with the jet grid. 
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However, at x=1.015, the source term model has a 
profile similar to the steady mass flow boundary 
condition without the jet grid farther away from the 
surface. As the flow field is examined further down 
field, at x=1.045, the source term model is almost 
identical to the steady mass flow boundary condition 
with out the jet grid. 
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Figure 20. u velocity distribution at x=0.985, v=0. 
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Figure 21. u velocity distribution at x=1.015. y=0. 
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Figure 22. u velocity distribution at x=l .030. y =0. 


X-1.045, y=0 



Figure 23. u velocity distribution at x= 1.045, y=0. 


Figures 24-26 show the u velocity contours at 
x=1.015 for the source term model and steady mass 
flow boundary condition cases. 
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Figure 24. u velocity contours of the source term 
model. 
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Figure 25. u velocity contours of the steady mass 
flow boundary condition without jet grid. 
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Figure 26. u velocity contours of the steady mass 
flow boundary condition with jet grid. 


The source term model is similar to the steady mass 
flow boundary condition with the jet grid close to the 
wall and to the steady mass flow boundary condition 
without the jet grid far from the wall. 
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Figure 27. v velocity distribution at x= 1.015 

Z=0.005. 
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Figure 28. w velocity distribution at x=1.015. 

Z=0.005. 


Figures 27 and 28 show the r and w velocity 
distribution at \=1.015 and z=0.005. respectively. 
The source term model does a fairly good job in 
modeling the shape of the v and w velocity 
distribution curves at the location specified. 
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Figure 29. Cross-stream vorticity. d . for source term 
model. 
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Figure 30. Cross-stream vorticity. d . for steady 
mass flow boundary condition without jet grid. 
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Figure 31. Cross-stream vorticity. d . for steady 
mass flow boundary condition with jet grid. 







The cross-stream vorticity, d , was calculated for all 
three cases. The vorticity is calculated using 
equation 5. 


-z -v 

Figures 29-31 show the cross-stream vorticity for all 
three cases at x= 1.015 (immediately after the jet). 
The source term model captures the vortical flow 
region fairly well, but the intensity is not as strong as 
the steady mass flow flow boundary condition with 
the jet grid. 

Conclusions 

In general, the source term model compared well 
with the steady mass flow boundary condition on 
both the 2D and 3D flat plates and is a viable 
alternative to gridding the jet. The source term 
model reproduces velocity profdes very similar to the 
steady mass flow boundary condition, which is used 
to model steady micro jets. Unlike the steady mass 
flow boundary condition, the source term does not 
need a grid for the jet to capture the angle of attack 
and incidence. This feature will allow an investigator 
to quickly investigate varying angles of attacks and 
incidence of a steady micro jet without having to 
generate a grid each time. 

The source term model is a viable tool for 
investigating steady micro jet installations on a 
surface. As large cases are investigated, the use of 
the model will significantly reduce the turnaround 
time from conception to solution by reducing grid 
generation and computational time. 
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